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On the basis of a proposed model of wave function collapse, we investigate spontaneous local- 
ization of a quantum state. The model is similar to the Ghirardi-Rimini- Weber model, while we 
postulate the localization functions to depend on the quantum state to suffer collapse. According 
to the model, dual dynamics in quantum mechanics, deterministic and stochastic time evolution, 
are algorithmically implemented in tandem. After discussing the physical implications of the model 
qualitatively, we present numerical results for one-dimensional systems by way of example. 
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I. INTRODUCTION 

There has been a significant increase of interest in the 
foundations of quantum mechanics (QM). This owes un- 
doubtedly to the technological progress achieved during 
the last decades in experimental investigation of how 
quantum systems behave not only statistically, but on an 
individual level. In fact, there are experimental results 
for microscopic systems which appear to be naturally 
explained by resorting to collapse of wave functionP, 
BlilllllElllilillSlil- Notwithstanding, wave 
function collapse has ever been the source of the prob- 
lem in understanding and interpreting QM, and there is 
still no definitive consensus on the "measurement" prob- 
lem. Accordingly, along with those which accept the 
collapse postulate, there are many other interpretations, 
such as de Broglie-Bohm causal interpretation 12 13 141 . 
statistical interpretationfisj . decoherenc efl^ lO . Il8lll9l 
I20I I2H . modal interpretation s l2p l24j . consistent 
(decoherent) histories 25, 26, 23, many- worlds/minds 
interpretations pM I291 130. 31, 331^ In this paper, we aim 
to investigate a realistic model of wave function collapse. 

First of all, it is totally unsatisfactory to connect wave 
function collapse with an act of measurement, because 
the concept of measurement is ill-defined [H, To 
describe wave function collapse from a realistic stand- 
point, we have to specify a well-defined dynamics of the 
stochastic process of collapse. In this regard, there are 
long-standing studies on stochastic nonlinear equations 
to realize collapse effectively. On the one hand, mod- 
els to cause rapid but continuous collapse in the sense 
of Brownian motion were studied by Bohm and Bubjs^, 
Pearlels^, Diosi'ST'l, and Gisin'ss'l. On the other hand, 
a model postulating discontinuous instantaneous collapse 
was elaborated in the pioneering work by Ghirardi, Rim- 
ini and Weber (GRW)!!! At present, the GRW 
model and the Pearle model have been jointly devel- 
oped to the continuous spontaneous localization (CSL) 
modelji^. Mathematically, CSL is based on a stochasti- 
cally modified Schrodinger equation. The current status 
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of studies on what are and have to be achieved in their 
models has been reviewed extensivelv |4^ liS li^. Ils j . On 
the basis of the idea of spontaneous collapse (SC), which 
may have a clue to the measurement problem, agood deal 
of related works are published indeed|4 ^ IIt I lia . liol l50l | . 
The task of theories based on and accounting for SC is to 
specify when, how and how often the collapse occurs, and 
to investigate the physical consequences. This is partic- 
ularly important because the effect of SC must be physi- 
cally relevant, in striking contrast to the other interpreta- 
tions. The predictions of collapse theories should more or 
less deviate from those of standard quantum mechanics 
(SQM). 

In this paper, we propose and investigate another such 
model of similar physical implications as the original 
GRW model. We also introduce two constants to char- 
acterize the dynamics of collapse. In the GRW model, a 
wave function is subjected to incessant spontaneous lo- 
calizations, or the wave function is multiplied by a local- 
ization function, e.g., the Gaussian function, of which the 
frequency and the localization width are postulated to be 
universal constants, A ~ 10~^^sec~^ and a ~ 10~^cm, re- 
spectively. The GRW localization mechanism is such that 
its frequency increases as the number N of constituents 
of a composite system increases, so that a macroscopic 
object comprising an Avogadro number of constituents 
collapses extremely rapidly, at a rate of NX ~ 10''sec~^. 
In contrast, we propose to postulate a constant collapse 
rate 70 = Tq^i while we model that the localization func- 
tions, particularly its length scale, are variable and de- 
pend on the wave function to suffer collapse. In effect, we 
introduce an energy scale Tq so that the stochastic col- 
lapse should occur to the effect to cost energy AE ~ Tq. 
Therefore, in our model, the macroscopic number A'' does 
not play any essential role, so that the effect of collapse 
may manifest itself even in microscopic systems |5l|. 

After introducing our collapse model in Section III Al 
we discuss the physical implications of the model in Sec- 
tion HTbI and in Section 1211 By way of illustration, in 
Section IIVI we apply the modified dynamics to obtain 
numerical results for one-dimensional systems. Although 
we are mainly interested in the individual dynamical evo- 
lution of a particular system, we briefly discuss the sta- 
tistical description of the system in terms of the density 
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matrix in Section We state our conclusion in Sec- 
tion lVll In Appendix^ we give an analytical treatment 
of a simple special case of the model to support general 
discussions in the main text. We discuss a generaliza- 
tion required to treat many particle systems in Appendix 
IbI The main purpose of this paper is to show that the 
proposed model provides us with a consistent picture of 
crossover between quantum and classical mechanics. 

II. MODEL 

A. Single particle problem 

We discuss a non-relativistic model of spontaneous col- 
lapse. Let us consider the normalized state ^'t(r) as a 
function of space r and time t, where the indices repre- 
senting internal quantum degrees of freedom like spins are 
suppressed for simplicity. We postulate that ^'t+At (r) af- 
ter the infinitesimal elapse of time At is determined from 
^t(r) stochastically as follows: (i) With the probability 

^t+Atir)= (^l-^HA?jMrl (1) 

while (ii) 

*t+At(r) = $„(r) = -Lp„(r)*t(r), (2) 

with the probability (1 — 7o)u'„Ai, where 

Wn = J dr|*t(r)|2F„(r)2. (3) 

In (PJ, the sum over the suppressed indices is implicit 
in the definition of the modulus square p(r) = |\['t(r)p, 
which thus represents the density in ordinary space. Un- 
like (i), the process (ii) is discontinuous by postulate, 
so that it cannot be put into a differential equation, 
though we can take the limit At ^ without affect- 
ing the physical content of the model. The branching 
structure of the proposed dynamics is schematically sum- 
marized in Fig. ^ According to the modified dynamics, 
the quantum state \E'(r) is subjected to collapse (ii) with 
the constant rate 70. The probability to keep following 
the deterministic evolution (i) for NAt = r is given by 

(1 — ^) exp (— 7ot) as — > 00. The collapse rate 
7o, or the time tq — 7(7 is the first parameter of our 
model. 

In |(5J, the effect of collapse is described by the real- 
valued continuous envelop functions Pn (r) , which we call 
the localization functions, following the GRW model. 
The factor ^yw^ in Q is to normalize the collapse out- 
comes, {^n\^n) — 1- The functions -Pn(r) have to satisfy 

^P„(r)2 = l, 0<P„(r)<l, (4) 

n 
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FIG. 1: Flow chart of the modified dynamics. From the 
initial state in the left, one (and only one) of the states in the 
right column is realized with the probability attached to the 
arrow. 



in order to meet the conservation of probability, 

J2w^ = l. (5) 

n 

To achieve spontaneous localization, the GRW model 
assumes P„ (r) to be the Gaussian functions with a fixed 
width centered at random positions. We postulate rather 
that a set of the real functions Pn(j) is determined op- 
timally depending on the state ^'((r) so as to maximize 

^' = 51^" /l*"l'log|*"l'dr, (6) 

n •' 

under the constraint 

ToAS = AE, (7) 

where 

A5 = -Tr(/51ogp) (8) 

and 

A£: = Tr(piJ) - (9) 

are the physical changes of entropy and energy due to 
collapse. They are defined customarily in terms of the 
density matrix p = pt+At of the collapse outcome. In 
real space representation, the matrix elements of p are 
given by 

(r|p|r') = ^w„$„(r)*<i>„(r'). (10) 

n 

The above is the procedure that we propose to fix 
the collapse outcomes $„ — resulting from an ar- 
bitrary state From our standpoint to regard wave 
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function collapse as a real process, ^ and (O represent 
the physical changes accompanied by the collapse, both 
of which are the important quantities to characterize the 
stochastic irreversible process. Thus {Tjl and Q respec- 
tively represent statistical conservation of 'free energy' 
F = E - ToS and the density. 



|vl/(r)p = ^|$„(r)|^ 



(11) 



Under these statistical constraints, we obtain the set 
of collapse outcomes by modifying to replace the 
squared amplitude p of the initial state = y/pe^" with 
Pn = P^p without affecting the phase a. Thus, in a 
sense, we envisage the collapse as a kind of fragmenta- 
tion of wave packet p — > {p„} in real space. 

In case where there happens to be non-equivalent sets 
of the localization functions {Pn} for given ^'t(r), e.g., by 
reason of symmetry, we require another postulate that a 
single set is chosen from among them at random with 
equal a priori probability. It would be clear that these 
postulates to formulate the probabilistic 'second law' (0) 
are adopted on the analogy drawn from Thermodynamics 
of the irreversible process of open systems immersed in 
a heat bath at temperature Tq. The constant Tq is the 
second parameter of our model. 

We adopt © as a proper expression to implement the 
spontaneous tendency of collapse to reduce the spatial 
extent of the quantum state, 

[Srf = ^i«„(($„|r.r|$„)-|($«|r|<i>„)P).(12) 

n 

Thus we take account of the special role of the position 
basis in terms of S". In effect, to minimize the variance 
((5r)^ means to maximize the second term in (I12|l . 



for the first term is independent of 



(13) 



w. 



($„|r •r|$„) = (*|r -rl^-). 



The expression in (|13|) is not easy to treat with formally. 
To serve our purpose of localization equally well, we shall 
rather adopt the expression ©. The positive sign in © 
is opposite to the conventional definition of entropy with 
a minus sign (cf. ©). This is because we intend to de- 
scribe the spontaneous tendency toward compact local- 
ization (instead of diffusion) of quantum states in terms 
of the physical quantity S' ~ log((^r)~^ to be maximized. 
For simplicity, for the same purpose, we might as well 
maximize a simpler expression like 



= ^w;„y dr|$„(r)|4. 



(14) 



Nevertheless, (0) has a physically preferable feature for 
the purpose of generalization (See Appendix IB|| . 



It is stressed that the physical implications of the pro- 
posed modified dynamics, which are discussed in what 
follows, are not affected qualitatively by whichever quan- 
tity one shall adopt to realize localization. The results 
of our primary concern are essentially the consequences 
of the very idea that stochastic collapse tends to reduce 
the spatial extent of wave packet as small as possible to 
the extent to cost energy AE ~ Tq by ||7J), or that the 
collapse outcomes are variably determined depending on 
the state ^' to be reduced. 



B. Quantum Limits 

By construction, the non-trivial effect of collapse be- 
comes negligible in the limit 70 — *■ 0. Consequently, in 
the time regime appreciably shorter than tq, the quantum 
results based on the deterministic evolution (i) remain in- 
tact. Hence, we may regard the short time regime i <C tq 
as a quantum regime. 

In the long run i ^ tq, however, the standard quantum 
predictions based on the strict validity of the unitary 
evolution (i) alone have to be modified more or less. One 
of the most notable consequences of the collapse model 
is the violation of energy conservation AE ^ ojsfll l4ll|. 
which is taken for granted in the above model through 
0. With the physical entropy AS" created upon collapse, 
the rate of energy production is given by TqAS/tq, which 
must be negligibly small in practice. 

For definiteness, let us discuss the system described by 
the Hamiltonian 



2m 



H 

for which we obtain 

ti2 



(15) 



AE ^ ^1 drp(r)5] (VP„(r))^ > 0, (16) 

n 

without approximation. It is noted that (|16|l depends on 
the potential V'(r) only implicitly through the 'probabil- 
ity' density 

P(r) = |*t(r)r 
In general, for an observable O, we obtain 
AO = ^zi;„($„|0|$„)-(*|6|*) 

n 

= /'drvI/,(r)*^P„(r)[d,P„(r)]_vl;,(r) (17) 

n 

Even when O and -Pn(r) commute, [0,P„(r)]_ = 0, the 
collapse ^ conserves the expectation value (O) only sta- 
tistically, AO = 0, so that we should generally expect 
(<I>„|d|$„) ^ (*|0|*) individually. For instance, the 
center of mass (r) will fluctuate as we follow the individ- 
ual behavior of a particular state evolving according to 
the modified dynamics (cf. Fig.O. 
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According to |7J and H16() . in the limit Tq 0, we 
conclude i-'ri(r) = constant if p{r) > for any r. There- 
fore, in this case, we recover not only the energy con- 
servation AE — 0, but also the quantum limit with no 
collapse effectively. This is because the collapse (0) with 
the constant P„ has no physical effect at all, for the col- 
lapse outcomes become physically all equivalent to the 
pre-coUapse state, l'„(r) = Pn'^ cx ^'(r). Hence, as well 
as 7o — !■ 0, the limit Tq ^ may also be called the quan- 
tum limit. Nevertheless, in contrast to the former limit, 
all the standard quantum results are not recovered in the 
latter. As discussed below in Section Iff! CI our modified 
dynamics eventually gives predictions against SQM even 
at To = 0. 



III. QUALITATIVE RESULTS 

A. Classical Limit 

According to the modified dynamics presented above, 
the spatial extent of the wave packet Vl/(r) is regularly 
truncated down to a finite size. Without detailed calcu- 
lation, we can make an order estimate of the length scale 
of the wave packet under spontaneous collapse. Owing 
to the constraint even in free space V{r) — 0, the 
energy scale Tq introduces the length 
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(18) 



which is nothing but the thermal de Broglie wavelength 
at "temperature" Tq. Accordingly, we conclude the finite 
wave packet width of order \6r\ ~ Ao cx m^^^^, which 
brings about our purpose to reproduce classical mechan- 
ics of the point mass r(t) = (r) in the macroscopic limit 
TO — !■ CX) where Ao — > 0, without introducing any external 
observer. The finite width 1)18(1 from is just as ex- 
pected by the thermodynamic analogy mentioned above. 
In effect, the degree of localization {5r)^ depends on the 
energy AE > required for the localization. If it were 
not for the constraint Q, or if we let Tg ^ oo, the hypo- 
thetical spontaneous collapse would reduce (Sr)"^ without 
limit, and the established results of quantum mechanics 
at a short length scale must be spoiled altogether. Thus, 
in the present model, the crossover length scale presumed 
between classical and quantum regimes is brought in by 
the non-trivial scale Tq in a controlled manner. 

The position (r) of the free particle fluctuates with the 
amplitude of the order of \Sr\ ~ Aq. In principle, the 
finite value of \5r\ should give rise to non-trivial quan- 
tum corrections to Newton's equation of classical me- 
chanics for the approximately well-defined position r{t) 
of the body. Nevertheless, it has to be remarked that 
the non-trivial effects due to \Sr\ ~ Aq can not nec- 
essarily be conspicuous in real situations, because they 
may be completely masked by environmental decoher- 
ence, the effects of which have been intensively discussed 



for decades '21', "55, '53^ . In effect , as we see below, the 
scale Tq of collapse should be substantially subtler than 
any perturbations of practical relevance. Still it would 
not be impractical to expect that the non-trivial predic- 
tions of the model, namely, the intrinsic decoherence of 
the constant rate 70, might survive the extrinsic distur- 
bances in some controlled situations. 



B. Born's rule 

Let us turn our attention to the microscopic system for 
which Ao defined in (|18|) is substantially larger than the 
typical length scale |(5r| ~ ct of the system. The latter is 
determined by the kinetic as well as potential terms in 
the Hamiltonian. 

To begin with, let us consider the wave packet whose 
linear dimension is of order a. We represent by A the 
typical length scale over which the localization functions 
Pn change appreciably. As we noted at the end of the 
last section, if p(r) > for any r, then we should obtain 
A ~ Aq cr. As a result, we conclude that the collapse 
has little effect on the quantum state, viz., Pn'i> cx 'J. For 
example, in the limit cr ^ Ao, the localized wave packet 
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(27rf72)3/2 



e-|r|V2. 



(19) 



is hardly affected by the postulated collapse. Hence the 
effect of (ii) may be schematically expressed as 



(jjair) 0^(r). 



(20) 



We conclude AE = TqAS ~ accordingly. 

Next, we consider the case in which the opposite limit 
A ^ cr may be realized, where the postulated collapse 
has a drastic effect. Since X]n(^^n('"))^ ((TCll takes a 
finite value ~ A"^ for r in the domain where P„ varies, 
we obtain 



AF~ 



(21) 



in which the factor q = /o(r)dr denotes the spatial 
integral over the 'boundary' of P„(r) where Pi{r)Pj{r) ^ 
for i ^ j (see (|A10|) '). Hence, by the constraint AE = 
TqAS — Tq, we obtain 



or 



A ~ Xoy/q, 



A ~ XqP, 



(22) 



(23) 



in terms of the density per length p defined by q = Xp. 
Therefore, A can be scaled down by the factor q or p, 
which may happen to be exponentially small, e.g., when 
p(r) is represented by a linear superposition of well sep- 
arated wave packets. 



For example, let us consider a typical case of an 'un- 
stable' linear superposition of the localized wave packet 

^-(0 =c+0<,(r-R) + c_0<,(r + R), (24) 

where |c+p + |c_p = 1. For *(r) with |R| > a, we 
find P±{v) to give P±^ ~ c±(t)a{v =F R-)- As a result of 
collapse, therefore, we obtain 

*^^0a(r-R), or -2^0,(r + R) (25) 

with the probabihties w+ = |c+p and W- = |c_p, re- 
spectively. The point is that the two states ^^(r + R-) 
and 0cr(r — R) have little spatial overlap, so that it hardly 
costs energy to reduce |^r| ~ 2(?ii+w_)^/^|R| » cr of 
the initial wave packet 1)24(1 down to \5t\ ~ a for (|25|l . 
which is 'stable' in the sense of (|20|l . Indeed, for def- 
initeness, owing to p cx e"'^' /'^ /cr, we obtain A <C (t 
for |R| 3> (T. Thus the localization functions are approx- 
imately represented by the step functions P_ (r) ~ 0(r-R) 
and P+{v) ~ 9{—Y ■ R), where 9{x) = 1 for x > and 
e{x) = for x < (cf. Figs. □ and IHl). This is essentially 
how Born's rule follows from our modified dynamics. 

The above argument is based on the emergence of the 
different scales A ^ (t <C Ao through |71l. As a function 
of the separation |R|, the qualitative change of behavior 
from H20|) to H25() occurs around |R| ~ cr, which is in- 
dependent of Tq. In other words, according to H22|l . the 
microscopic scale A for P„ may be derived as the prod- 
uct of Ao, which must be larger than any other length 
scales of a microscopic quantum system, times the expo- 
nentially small factor y^. Although the former depends 
on To , it is essentially the latter factor q that provides us 
with the smallest scale A of P„. Thus the qualitative fea- 
tures discussed above do not depend sensitively on the 
precise value of Tp, but are mainly determined by the 
density distribution of the state ^'(r). 

C. Predictions against SQM 

The modified dynamics describes the 'measurement' 
involving wave function collapse as a stochastic process 
to minimize \5v\ of the wave function It serves our 
purpose that those classically inadmissible superposition 
states of well separated wave packets like (I24|l are desta- 
bilized to decay with the finite lifetime tq. In princi- 
ple, the non-trivial physical effects of the model should 
be manifested themselves in the interference experiments 
which are devised to distinguish the pure state 5* from 
the statistical mixture of the collapse outcomes $„. For 
instance, we are led to an unconventional prediction of 
the destruction of the quantum interference pattern, as 
schematically presented in Fig. |2] Since such an intrin- 
sic decoherence is concluded independently of the actual 
distance between the separated wave packets in super- 
position, we must also conclude seemingly unfavorable 
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FIG. 2: According to the modified dynamics, the interfer- 
ence fringe pattern due to split and recombined wave packets 
is predicted to be washed out, if the distinct wave packets in a 
superposition had been spatially separated further than their 
own widths for a duration longer than the collapse time ro 
before the recombination, e.g., if the length I of the impene- 
trable partition slab separating the two paths is made larger 
than a coherence length ~ «ro, where v is the velocity of the 
incoming wave packet. 



consequences for the microscopic single-particle states 
with nodes along which the density vanishes, p(r) — 0. 
In contrast to (|20|) for (|19|l . e.g., the wave function 
07r (r) oc x(j)cr (f) for cr ^ Ao cannot but be halved to result 
in the effectively nodeless collapse outcomes P±(t)n] 

as one obtains A ~ ct'^/Aq ^ in this case. In short, 
the modified dynamics predicts that wave functions with 
nodes generally have to be nodeless eventually. It is clear 
that such drastic unconventional collapse processes must 
not intervene the conventional Schrodinger time evolu- 
tion too frequently. Thus we have to require the collapse 
time To to be substantially slower than any practical co- 
herence time scales for the quantum interference of spa- 
tially separated wave packets. 

D. Numerical estimate of the model parameters 

In practice, the numerical value for To must be prop- 
erly chosen such that the characteristic length Aq for mi- 
croscopic particles like electron and neutron must be too 
large while the intrinsic wave packet widths Ao for macro- 
scopic bodies, which may be effectively regarded as ma- 
terial points, are too small to be unexpectedly detected 
experimentally. To put it concretely, in order to guaran- 
tee Ao > 10m for electron while Ao < l/.tm for a tiny par- 
ticle of a nanogram, we estimate 10~^^ J < To < 10"'^^ J. 
Hence there is a relatively wide range open for To, owing 
to the above note that the model is relatively insensitive 
to To. 
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On the other hand, we must assume the collapse time 
To to be larger than a typical coherence time scale of mi- 
croscopic quantum phenomena, as mentioned in the last 
subsection. Nonetheless, at the same time, tq must be 
shorter than a typical scale of macroscopic classical phe- 
nomena in order to warrant definite outcomes for classical 
bodies approximately at any time, i.e., to prevent such 
an embarrassing superposition state of 'classical' states 
with 'distinct' spatial configurations from developing and 
lasting for a long while. Therefore one may make a rough 
estimate like 10~^sec ^ tq ^ 10~"'^sec, so we find that 
the parameter range allowed for tq is more restricted than 
Tq. In fact, we consider that whether we could assume tq 
properly would be the crucial point for the present model 
to be viable physically. In what follows, we discuss the 
implications of the model assuming that experiments do 
not exclude room for the model parameters. 

Without being affected by the collapse, the wave 
packet evolves according to ([[J for the duration of or- 
der Tq. In the meantime, the width (i5r)^ ~ Ag for a free 
particle grows up to Aq (l -I- i^^)"^)), so that the second 
length scale must be introduced if To'^'o ^ 47r?i. However, 
with the above estimate of Tq taken for granted, this limit 
must be rejected physically, because we hardly accept tq 
as long as h/To > 10^ second on physical grounds. Hence 
it is valid to use the single length scale 1)18(1 in free space, 
since we can consistently assume 

TqTo < Anh (26) 

without invalidating the model. 



IV. MOTION IN ONE DIMENSION 
A. Model 



ization functions as follows (cf. Fig.|Hl), 

{1, X — xq < — 2A 

cos^(a;-a;o-f 2A), |a;-xo|<2A (28) 

0, 2X<x-xo 

and 

( 0, x-xq<~2X 

P^o(x-) = < cos^(x- xo - 2A), |a;-a;o|<2A (29) 

I 1, 2X < X — xq 

with which we obtain 

The entropy change is approximately given by 

AS -WRlogWR - WLlogWL, (31) 

where 

/oo 
{Pl°{x)fp{x)dx 
-oo 

and WR = 1 — WL are the probabilities to result in P^^vj/ 
and P^°'i', respectively. 

In order to fix xq and A, we have to evaluate S' in ® 
for the given density p{x) = |^'(a;)p at every moment. 
Nevertheless, to reduce the numerical task, we assume to 
set A = Ao/8 arbitrarily, and choose xq so as to minimize 
AF = AE - TqAS with (ETJ in anticipation that the 
qualitative features shown below would not be severely 
modified by the assumption. In effect, we are interested 
in the interplay of the stochastic and deterministic evo- 
lution. In what follows, we set h = 1, while assuming 
To = 1 and tq = 1 despite (|^ . 



To give shape to the model, we investigate motion of a 
particle in one dimension. The unmodified Schrodinger 
equation may be written in a dimcnsionless form. 



d 



(27) 



where t — TqI/K and x = x/Xq. As discussed above, for 
macroscopic classical bodies, the potential V(x) should 
be slowly varying, |VI^| <C \V\. In practice, the time 
scale To /hot (|?7| should be substantially slower than tq, 
as noted in (|26|l . Nevertheless, in this section, we assume 
TqTo = h and t — I/tq, for simplicity, to illustrate how the 
probabilistic and deterministic time evolutions compete 
with each other. 

It is generally difficult to evaluate the localization func- 
tions analytically for an arbitrarily given state (cf. Ap- 
pendix^, so that we discuss only the simple case that 
the state '^{x) is divided by collapse into the left and 
right halves at a position xq. To describe the collapse, 
let us make use of the two trial functions for the local- 



B. Results 

To begin with, we give results for a free particle, 
V{x) — 0, in Fig. 131 Here and below, the results for the 
density p{x, t) with and without collapse are compared 
with each other to make the non-trivial collapse effect 
clear. As shown by contours on the base {x-t plane) of 
Fig. O the collapsing wave packet traces a zig-zag world 
line. To a greater or lesser extent, the discontinuity is 
a general and non-trivial feature to be revealed in our 
results. This is made conspicuous when we should oth- 
erwise, without collapse, get involved in a controversial 
superposition state. As a typical case, results for a mov- 
ing particle impinging on a potential barrier at x ^ are 
given in Fig. ^ 

Results for the system in a double well potential are 
shown in Fig. [SJ for which we varied xq for fixed A — 
Ao/16, as the constraint allows the width A to be smaller 
in this case. We find that even A = Ao/32 does not 
change the results qualitatively, for the wave function 
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FIG. 3: Time evolution of the wave packet density p(x,t) of a free particle with 2m = 1, starting from the Gaussian 
5'(a;,0) oc e"^^''*' with S = Ao/3 (cf. Ao = 3.5). Without (left) and with (right) collapse. Two stochastic collapses are 
recognized. In our model, the intermittent collapses keep the wave packet from diffusing, while setting it in Brownian motion. 




FIG. 4: Quantum tunneling through a potential barrier at a; ~ 0. The density p(x,t) of a particle with m = 1 in motion, 
evolving from *(x,0) oc e-(^/'')''+'fo^o with 5 = Ao/2, xo = -8 and po = 8. The potential is V{x) ~ O.OSe"''^' . Without (left) 
and with (right) collapse. Two collapses before and after refiection(transmition) are recognized. Tunneling failed, and the wave 
packet is reflected. Aside from a digression over the barrier, the wave packet traces an approximately well-deflned trajectory. 



tends to collapse around x ~ 0. The suppression of a 
coherent quantum motion by frequent "measurements" is 
generally known as the quantum Zeno effect [H^ Issl l56l| . 
This is demonstrated in Fig.[S| We are also interested in a 
related problem of quantum decay from a potential well. 
To realize such, we solved a simple model as indicated 
in Fig. It is noted that the unstable state trapped 
in the potential barrier looks stationary until it decays 
spontaneously. 

Lastly, to mimic situations encountered in position 
measurement, wave packet reduction induced by a time- 
dependent potential is displayed in Fig. [7| which shows 
that gradual increase of the external potential may cause 
a sudden change of state. To obtain Fig. [7| we varied A 
for xq = 0. This behavior corresponds to the transition 
from (|20|) to (|25|l . In general, we associate the "measure- 



ment" of a quantum particle with such an unordinary 
leap caused by a specially prepared external potential. 



These results demonstrate that the effect of collapse 
becomes conspicuous when the linear width 5 of the wave 
packet exceeds the characteristic length scale Aq, and 
that the collapse hardly affects the wave packet when 
(5 <C Aq. Therefore, as discussed in the previous sections, 
we obtain a picture in that quantum (wave) mechanics 
obtains in the one limit Aq ^ oo, while classical (parti- 
cle) mechanics (with a definite path) emerges in the other 
limit Aq — > 0. 
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FIG. 5: Quantum oscillation in a double well potential (m = 2). Without (left) and with (right) collapse. There are potential 
walls at a; = ±1, and a square potential of the width and the height of 0.1 at x = 0. The sinusoidal continuous oscillation 
apparent in the left figure is suppressed by collapse. On the time scale of the figure, individual collapse are not conspicuous. 
In the right figure, we observe that the particle is found to be either in j: > or a; < 0, instead of being in their superposition. 
The particle position changes discontinuously from time to time at random. In a conventional phrase, the particle behaves as 
if it is "measured" with the frequency 70. 





FIG. 6: Quantum decay of a particle m = 2 trapped in a potential well. The potential well —1 < 2; < is that of Fig.^with 
the wall at a; = 1 removed. Without (left) and with (right) collapse. The wave packet is literally trapped until it goes off at a 
definite, but unpredictable, instant. 



V. DISCUSSION 

According to the model, wave function reduction 
causes the density matrix of the pure state vE'(r) to evolve 
as 

p(r,r') ^'(r)^'*(r') — > Z\(r, r')*(r)**(r'), (32) 
where 

zi(r,r')^5]P„(r)P„*(r') (33) 

n 

equals 1 for r = r' because of (0J, and vanishes for 
|r — r'l ^ A. Hence, by postulate, the diagonal ele- 
ments p(r, r) = p(r) are invariant upon collapse. As 



the functions P„ are determined by p(r), so is Z\(r, r'). 
Therefore, H32I) generally represents a non-linear trans- 
formation. Nevertheless, in many practical cases, the lo- 
calization functions may be determined by the external 
potential V{r), as in Figs. |S1 and [T] In fact, this would 
be the case when we can neglect wave function collapse 
except in measurement-related situations. In such cases, 
we obtain a linear equation, referring to Fig. ^ 

pt+dt = {I - lodt) (^pt + ^[pt,H]^dt^ +7odtZ\(r,r')pt, 
or 

^ = i[p,i/]_-7o(l-^)p. (34) 
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FIG. 7: Quantum "measurement" caused by a time- dependent potential. A partition wall is adiabatically inserted at the 
center of a bound state (2m = 1) in the potential well —1 < a; < 1. Without (left) and with (right) collapse taken account of. 
The potential is the same as in Fig.|K| except that the height of the central square potential is increased as ^height = O.li. The 
wave packet is finally found in either of the left or right partition. 



The equation of motion (|34|) has the same form as that 
investigated by Joos and Zeh to assess the effect of deco- 
herence through interaction with the environment 18]. 

At the statistical level of description, the collapse 1)32(1 
does not affect the density matrix p{r,r') unless it has 
the longer ranged correlation than the correlation length 
A of H33|) . For example, for free particles in thermal equi- 
librium at temperature T, we have 



p(r,r') ^ A^^/'exp (- '''\/'' ) , (35) 



where At is the thermal de Broglie wavelength at T. In 
order for H32I) to affect the statistical results based on 
(|35|l . we should assume At ^ Ao, or the physical tem- 
perature must be extremely low, T <C To- Consequently, 
in the statistical ensemble, the phenomenological effects 
of collapse apparent at the individual level can be com- 
pletely masked by environmental decoherence [20l l2ll Is^ 
l53. 57]. In other words, the discontinuity due to collapse 
as indicated in Figs. El and [3 would be averaged out at 
the statistical level. 

Last but not least, the theory discussed throughout the 
paper is essentially non-relativistic. This is evident from 
a special role we assigned to time t. It is still noteworthy 
that the model formulated with p{r) is Galilei invariant. 
The conflict of the instantaneous collapse with special 
relativity has been elucidated [sl]- In this regard, we only 
note a consistent if naive attitude to fit with the present 
model, that is, to assume a special frame of reference, as 
noted by Eberhard'5^ and BcU'gS'], in which the collapse 
dynamics applies specifically 61] . 



VI. CONCLUSION 



In this paper, we addressed ourselves to the problem 
of how the classical world as we experience it emerges 
from the underlying laws of quantum mechanics. We 
investigated the motion of a quantum particle on the ba- 
sis of the proposed modified dynamics of wave function 
collapse. In effect, the proposal may look complicated 
in practice, owing to the task of solving for the local- 
ization functions, but the model is simple in principle, 
as summarized in Fig. ^ We describe the collapse dy- 
namics in terms of physical quantities such as energy, 
entropy, and density in real space. The model contains 
two material- independent constants tq — 'Jq^ and Tq. It 
is suggested on the analogy of the second law of ther- 
modynamics in order to destabilize those linear super- 
positions of spatially separated wave packets, which are 
hardly interpretable from an ordinary, classical, realistic 
standpoint. On the basis of the proposal, the numeri- 
cal results are obtained as displayed in Figs. I3I7I The 
model reproduces classical mechanics as a limit, yet it 
does not use the classical limit for its own formulation. 
According to the model, we obtain a unified mechanical 
picture in that the classical picture emerges on a long 
and large scale, t ^ tq and a; ^ A, while the quantum 
picture applies in the opposite limits, t <^ tq and x <C A, 
where A denotes the natural linear dimension of the wave 
packet determined according to the model. In fact, the 
length scale A of collapse depends on circumstances, and 
it can be made arbitrarily small. This is particularly the 
case for a microscopic particle subjected to measurement- 
related situations, where one can devise an apparatus to 
make A as small as one likes, in principle. 
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APPENDIX A: VARIATIONAL PRINCIPLE 

In principle, it is possible to write down the equations 
for P„(r), that is, the equations to determine the collapse 
outcome states ^^(r). However, it is generally difficult 
to solve them for an arbitrary state ^'(r) analytically. In 
this appendix, we investigate the simplest tractable case 
where a wave packet is justly halved with equal proba- 
bility by collapse. 

We obtain the localization functions Pi and P2 under 
the condition that the collapse probabilities to the two 
outcomes are equal, viz., wi ~ W2 ^ 1/2. To this end, 
the auxiliary function 9(r), ranging from to 7r/2, is 
introduced by Pi = cos 6 and P2 = sin 9 to take account 
of the constraint P^ + P2 — 1. In general, the collapse 
outcomes Pn^* are not orthogonal, or the overlap integral 



(4-1^1 P2I*) = ^ J drp(r)sin26l(r) (Al) 

does not vanish. Hence the entropy production by the 
collapse is given by 



A^ = -(i+p)log(i+p)-(l-p)log(i-p 



which takes the maximum value log 2 for p = 0, while the 
minimum is AS" = for p = 1/2. 

In order to achieve the localization of wave function, 
to maximize in 114() is more easily implemented than 
to work with ((5r)^, (|12(l . Hence, to start with, let us 
maximize 



where A = ^JW^yfj^mp. 

The equation of a pendulum (jASp can be solved ana- 
lytically. The solution with a proper boundary condition 
is given by 

cos W = -sn(x/A, 1) = - tanh(a;/A), (A6) 

where sn(u, fc) is the Jacobian elliptic function. The ori- 
gin of the X axis has to be consistently chosen to give 
wi = W2 = 1/2, i.e., by 

y tanh(x/A)pdr = 0. (A7) 

The solution ljA6|l for A = 1 is shown in Fig. |S1 with the 
solid curve. Note that the curve essentially represents 
the localization function P^ixY ■, because 

„ , ,2 1± 00826* 

PiA^r = — ^ — • 

To estimate the energy cost (|A3|) to fix A, we can make 
use of 'energy conservation' of the pendulum motion 
"[1, 



8A^(^) +V{e)^V{0), 



where 



V(e) = - cos 46*. 
4 

From the last equations, we obtain 



(A8) 



(A9) 



P^ = '2J (sm^ + cos^ 6l)p(r)2dr (A2) 
under (TJ), for which 

AE ^ ^ f {Vef p{r)dr. (A3) 
2m J 

In terms of a Lagrange multiplier A', the equation for 9 
is obtained by the variation principle, 

A(^+A'(ToA5-Ai?)) = 0, 
from which we obtain the equation 

A' ( —V(pVe) -To^ log 008 26*^ = 2p'^ sin 46*. 
\m 2 1 — p J 

(A4) 

In particular, we are interested in the case p ~ and 
|Vp/p| <C |V6'|, where the model predicts the non-trivial 
result that the collapse outcomes Pi 4" and P2* are dis- 
tinctly different from each other, ($i|$2) — 0. To discuss 
this case further in depth, regarding as a function of x 
for simplicity, we obtain the one-dimensional equation 



4A2^ = sin 40, 



(A5) 



4AM^) ^PiixfP^ixf 



hence. 



where 



AE = 



8mA2' 



(AlO) 



q = j P^{x)Pi{x)pdr ~ J e^^'/^Vdr. (All) 

With respect to the last approximation, we should note 
J Pi{x)'^P2{xfdx = 2A, while /e-^'/^'da; = V^A, and 
Pi{x)^P2{x)^ ~ e-^'/^' ~ 1 - (x/A)2 for \x\ < A. Simi- 
larly, we may write (|A1|) as 



p : 



pdr. 



(A12) 



with which p(v) — |4'(r)p must meet p ^ 1 as well 
as (|A7ll for consistency. Given such p(r), the parame- 
ter A for Pn(r) is determined by AE ~ Tolog2, as in 
(|22Jl. In consequence, we obtain the collapse outcomes 
<f>„ = \/2Pn^- Note in particular that P„ becomes the 
step function in the limit A — > 0, e.g., as a result of the 
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FIG. 8: The localization function 2P2{xf - 1 = - cos26i is 
shown as a function of x. The solid and dashed curve are 
the solutions of l|A5|l and l|A13|l for A = 1, respectively. The 
dotted curve is a sine curve, which is added for reference. 



low density at x = 0, p{x = 0) ~ 0. Then we recover 
the Born's-rule probabilities Wn — |(<i>„|^)P together 
with the non-overlapping collapse outcomes owing 
to ~ Pn- As is clear from the qualitative discussion 
in Section flll Bl this feature is generic and not specific to 
the assumption wi = W2. 

In a similar manner, we can fix P„ by maximizing © 
instead of p^, HA2|I . To obtain <!>„ from our task is to 
maximize 

S' + (3iToAS- AE), 

where E = TripH), S = -Tr(/51ogp), and (3 is the 
Lagrange multiplier to keep the postulated constraint 
A(E — TqS) = 0. In the same approximations as above, 
in place of (jA5|) . we obtain 

— sin 29 log cot 6*, ( A13) 

dx^ 

where 4A^ — H^fi/m. A numerical solution cos26' for 
A = 1 is shown in Fig. |H1 with the dashed curve. We 
find that the optimal shape of the localization function 
Pn{x) is not severely modified. For reference, a sinusoidal 
fimction sin7ra;/4 for \x\ < 2 is also shown in Fig.|Slwith 
the dotted curve. 

It is noted that the solution of IjAlSp gives a finite 
value for \x{9 = 7r/2) — x{6 — 0)|, which is the over- 
lap width of Pi{x)P2{x). This is because the 'potential' 
V{6) = (logsin(26') + cos(26l)logcot6')/2 for (|XT3^ makes 
the integral 

Jo ^vio) ~ Vie) 

convergent, while ljA9p makes this expression divergent. 
Thus we obtain compact support of $„, which is a favor- 
able feature of the model. Moreover, from (jA8|) . we ob- 
tain |d6'/dx| — > as 6* ^ 0, 7r/2 for the solution of HA13|I . 



by which follows the continuity of $„(a;) as well as their 
derivatives as a function of x. Furthermore, as discussed 
below, the expression © is generalized adequately to de- 
scribe independent collapse of many particle systems. 



APPENDIX B: GENERALIZATION TO 
MANY-PARTICLE SYSTEMS 

In the main text, we focused on single-particle prob- 
lems of the modified dynamics. To describe collapse 
of a many particle system, we have to generalize the 
localization functions to P„(ri,r2,-- - ,r7v), which are 
to be multiplied with the N particle wave function 
^(ri,r2,--- ,rff). Accordingly, the integral / dr in Q 
is generalized to / dridr2 ■ ■ - drN over 3A^ dimensional 
space for the N particles. 

For N distinguishable particles, one should maximize 

S' = Y,Wn f |$„|Mog|$„pdridr2---drAr, (Bl) 

instead of In particular, when the N particles are 
independent from each other, the wave function of the 
whole system is represented by the product of wave func- 
tions, 

N 

5'(ri,r2, • • • ,rAr) = J| 

for which one can regard formally that each of the states 
■01 (fi) collapses independently, 

since (|Blll is decomposed into a sum of N independent 
terms, 

N 

S' - EE^". / I0«j'log|0„j2dr„ 
1=1 fii •' 

owing to Wn — Y\i Wm ■ This is a feature to be expected 
physically. 

For indistinguishable particles, the symmetry of the 
wave function ^' must be respected, so that the local- 
ization functions P„(ri,r2,- • • ,rAr) must be symmetric 
under the exchange of the arguments. For definiteness, 
they are approximately given in terms of P„ (r) for a sin- 
gle particle state as follows. The simplest assumption is 
to regard them simply as functions of the center-of-mass 
coordinate R = Vi/N, 

P„(ri,r2, • • • ,rAr) = P„(R). 

Otherwise, they may be represented as symmetrized 
products of P„ (r) , 

P„(ri, r2, • • • , rAr)2 = ^ P„, (ri)2p„^ (r2)2 • • • P„„ (rjvf, 
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where n represents a set of integers (rii,n2, ....,nj^) and 
the sum in the right-hand side is taken over all different 
permutations of them, so that 

n i—1 \ Ui I 

For identical particles, the following expression may 
equally serve our purpose, 

S' = Y^Wn\ Pn{r) logp„(r)dr, (B2) 



where /o„(r) = ($ri|'^r|'^'ri) denotes the particle density 
in terms of the number operator hy.. In any case, the 
identical particles must evolve simultaneously as 4* ^ 
Pn'^ / y/wZ, as they are not strictly independent from each 
other. 

To treat the case where the particle number N is vari- 
able, we have to generalize the collapse prescription fur- 
thermore to realize the localization on the number basis 
as well. This is not difficult in principle, but goes beyond 
the scope of the paper. 
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